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We derive an expression for a density operator estimated via Bayesian quantum inference in 
the limit of an infinite number of measurements. This expression is derived under the assumption 
that the reconstructed system is in a pura state. In this case the estimation corresponds to an 
averaging over a generalized microcanonical ensemble of pure states satisfying a set of constraints 
imposed by the measured mean vàlues of the observables under consideration. We show that via the 
"purification" ansatz, statistical mixtures can also be consistently reconstructed via the quantum 
Bayesian inference scheme. In this case the estimation corresponds to averaging over the generalized 
canonical ensemble of states satisfying the given constraints, and the reconstructed density operator 
maximizes the vonNeumann entropy (i.e., this density operator is equal to the generalized canonical 
^ ■ density operator which follows from the Jaynes principle of maximum entropy). We study in detail 

the reconstruction of the spin-1/2 density operator and discuss the logical connection between the 
three reconstruction schemes, i.e., (1) quantum Bayesian inference, (2) reconstruction via the Jaynes 
principle of maximum entropy, and (3) discrete quantum tomography. 
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I. INTRODUCTION 



The essence of the problem of state determination (= state reconstruction) lies in an a posteriori estimation of 
a density operator (= corresponding quasiprobability density distribution) of a quantum-mechanical (microscopic) 
system based on data obtained with the help of a macroscopic measurement apparatus ||ï[. The quality of the 
reconstruction depends on the quality of the measured data and the efficiency of the reconstruction procedure with 
the help of which the inversion-data analysis is performed. In particular, we can specify three different situations: 
Firstly, when all system observables are precisely measured; Secondly, when just part of the system observables is 
precisely measured; Finally, when measurement does not provide us with information enough to specify exact mean 
' vàlues (or probability distributions) of observables under consideration. 



A. Complete observation level 

Providing all system observables (i.e., the quòrum ^) have been precisely measured, then the density operator of a 
quantum-mechanical system can be completely reconstructed (i.e., the density operator can be uniquely determined 
based on the available data). In principle, we can consider two different schemes for reconstruction of the density 
operator (or, equivalently, the Wigner function) of the given quantum-mechanical system. The difference between 
these two schemes is based on the way in which information about the quantum-mechanical system is obtained. The 
first type of measurement is such that on each element of the ensemble of the measured states only a single observable 
is measured. In the second type of measurement a simultaneous measurement of conjugated observables is assumed. 
We note that in both cases we assume an ideal, i.e., unit-efficiency, measurements. 



*We dedicate this paper to the sixtieth birthday of Professor Jan Pefina. 
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1. Quantum tomography 



When the single-observable measurement is performed, a distribution W\^^{A) for a particular observable A in the 
State is obtained in an unbiased way i.e., — where \^a) are eigenstates of the observable À 

such that J2a 1=1- Here a question arises: What is the smallest nuniber of distributions W\^) [A] required to 

determine the state uniquely? If we consider the reconstruction of the state of a harmònic oscillator, then this question 
is directly related to the so-called Pauli problem Q of the reconstruction of the wave-function from distributions 
W\^,){q) and W^|if)(p) for the position and momcntuni of the state As shown by Gale, Guth and Trammel Q for 
example, the knowledge of W\-^\i{q) and VFj^^(p) is not in general sufficient for a complete reconstruction of the wave 
(or, equivalently, the Wigner) function. In contrast, one can consider an infinite set of distributions W\^){xg) of the 
rotated quadrature xe = qcos0 +psm9. Each distribution W\^i/^{xg) can be obtained in a measurement of a single 
observable íe, in which case a detector (filter) is prepared in an eigenstate \xe) of this observable. It has been shown 
by Vogel and Risken |^ that from an infinite set (in the case of the harmònic oscillator) of the measured distributions 
W\^í,){xe) for all vàlues of 9 such that [0 < 6* < tt], the Wigner function can be reconstructed uniquely via the inverse 
Radon transformation. In other words knowledge of the set of distributions W\iii'^{xg) is equivalent to knowledge of 
the Wigner function (or, equivalently, the density operator). This scheme for reconstruction of the Wigner function 
(the so called optical homodyne tomography) has recently been realized experimentally by Raymer and his coworkers 
. In these experiments Wigner functions of a coherent state and a squeezed vacuum state have been reconstructed 
from tomographic data. Quantum-state tomography can be applied not only to optical fields (harmònic oscillators) 
but for reconstruction of other physical systems, such as atòmic waves (see recent work by Janicke and Wilkens 
Leonhardt |^ has recently developed a theory of quantum tomography of discrete Wigner functions describing states 
of quantum systems with finite-dimensional Hilbert spaces (i.e., angular momentum or spin). 



2. Filtenng with quantum rulers 

In the case of the simultaneous measurement of two non-commuting observables (let us say q and p) it is not 
possible to construct a joint eigenstate of these two operators, and therefore it is inevitable that the simultaneous 
measurement of two non-commuting observables introduces additional noise (of quantum origin) into measured data. 
This noise is associated with Heisenberg's uncertainty relation and it results in a specific "smoothing" (equivalent 
to a reduction of resolution) of the original Wigner function of the system under consideration (see Refs. and 
[pl|). To describe a process of simultaneous measurement of two non-commuting observables, Wódkiewicz |ï^ | has 
proposed a formalism based on an operational probability density distribution which explicitly takes into account the 
action of the measurement device modelled as a "filter' (quantum ruler). A particular choice of the state of the ruler 
samples a specific type of accessible Information concerning the system, i.e., Information about the system is biased 
by the filtering process. The quantum- mechanical noise induced by filtering formally results in a smoothing of the 
original Wigner function of the measured state ||ï^,|rï| , so that the operational probability density distribution can be 
expressed as a convolution of the original Wigner function and the Wigner function of the filter state. In particular, if 
the filter is considered to be in its vacuum state then the corresponding operational probability density distributions 
is equal to the Husimi (Q) function The Q function of optical fields has been experimentally measured using 

such an approach by Walker and CarroU p3[ |. The direct experimental measurement of the operational probability 
density distribution with the filter in an arbitrary state is feasible in an 8-port experimental setup of the type used 
by Noh, Fougéres and Mandel |Q (for more detalls, see the recent book by Pefina and coworkers p^). 

As a consequence of a simultaneous measurement of non-commuting observables the measured distributions are 
fuzzy (i.e., they are equal to smoothed Wigner functions). Nevertheless, if the detectors used in the experiment have 
a unit efficiency (in the case of an ideal measurement), the noise induced by quantum filtering can be "separated" 
from the measured data and the density operator (Wigner function) of the measured system can be "extracted" from 
the operational probability density distribution. In particular, the Wigner function can be uniquely reconstructed 
from the Q function. This extraction procedure is technically quite involved and it suffers significantly if additional 
stochastic noise due to imperfect measurement is present in the data. 

We note that propensities, and in particular Q-functions, can be also associated with discrete phase space and 
they can in principle be measured directly p6| . These discrete probability distributions contain complete Information 
about density operators of measured systems. 
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B. Reduced observation levels and MaxEnt principle 



As we have already indicated it is now well understood that density operators (or Wigner functions) can, in principle, 
be uniquely reconstructed using either the single observable measurements (optical homodyne tomography) or the 
simultaneous measurement of two non-commuting observables. The completely reconstructed density operator (or, 
equivalently, the Wigner function) contains information about all independent moments of the system operators. For 
example, in the case of the quantum harmònic oscillator, the knowledge of the Wigner function is equivalent to the 
knowledge of all moments ((à^')™à") of the creation (à^) and annihilation (ò) operators. 

In many cases it turns out that the state of a harmònic oscillator is characterized by an infinite numbcr of inde- 
pendent moments ((à^)'"a") (for all m and n). Analogously, the state of a quantum system in a finite-dimensional 
Hilbert space can be characterized by a very large number of independent parameters. A complete measurement of 
these moments can take an infinite time to perform. This means that even though the Wigner function can in principle 
be reconstructed the coUection of a complete set of experimental data points is (in principle) a never ending process. 
In addition the data processing and numerical reconstruction of the Wigner function are time consuming. Therefore 
experimental realization of the reconstruction of the density operators (Wigner functions) for many systems can be 
problemàtic. 

In practice, it is possible to perform a measurement of just a finitc numbcr of independent moments of the system 
operators which means that only a subset Gi, {v = 1, 2, n) of observables from the quòrum (this subsct constitutes 
the so-called observation level ||ï^) is measured. In this case, when the complete information about the system is not 
available, one needs an additional criterion which would help to reconstruct (or estimate) the density operator uniquely. 
Provided mean vàlues of all observables on the given observation level are measured precisely, then the density operator 
(the Wigner function) of the system under consideration can be reconstructed with the help of the Jaynes principle of 
maximum entropy (the so called MaxEnt principle) |r^]. The reconstructed density operator fulfills several conditions. 
Firstly, its trace has to be equal to unity (i.e., Trp = 1). Secondly, Tr(/5G^) = Gi, {v = 1,2, ...,n) which means that 
the reconstructed density operator provides us with the measured mean vàlues of those observables which constitute 
the given observation level. Obviously, a large number of density operators can fulfill these two constraints. So one 
needs an additional criterion which would uniquely specify the generalized canonical density operator. According to 
Jaynes this operator has to be that one with the largest value of the vonNeumann entropy S — — Tr{/51n/5}. This 
additional condition means that the MaxEnt principle is the most conservative assignment in the sense that it does 
not permit one to draw any conclusion not warranted by the experimental data . 

The MaxEnt principle provides us with a very efíicient prescription how to reconstruct density operators of quantum- 
mechanical systems providing mean vàlues of a given set of observables are known. It works perfectly well for systems 
with (semi)infinite Hilbert spaces (such as quantum-mechanical harmònic oscillator) as well as for systems with 
finite-dimensional Hilbert spaces (such as spin systems). If the observation level is composed of the quòrum of the 
observables (i.e. this is a complete observation level), then the MaxEnt principle represents an alternative to quantum 
tomography, i.e. both schemes are equally suitable for the analysis of the tomographic data (for detalls see |Q). To 
be specific, the observation level in this case is composed of all projectors associated with probability distributions of 
rotated quadratures. The power of the MaxEnt principle can be appreciated in analysis of incomplete tomographic 
data (equivalent to a reconstruction of the Wigner function in the discrete phase space). In particular, Wiedemann pC| ] 
has performed a numerical reconstruction of the density operator (Wigner function) from incomplete tomographic 
data based on the MaxEnt principle as discussed by Buzek et al. |jï^. Wiedemann has shown that in particular 
cases MaxEnt reconstruction from incomplete tomographic data can be several orders better compared to a Standard 
tomographic inversion. These results can be interpreted as suggesting that the MaxEnt principle is the conceptual 
basis behind the tomographic reconstruction (irrespective whether in continuous or discrete phase spaces). 

C. Incomplete measurement and Bayesian inference 

It has to be stressed that the Jaynes principle of maximum entropy can consistently be applied only when exact 
mean vàlues of the measured observables are available. This condition implicitly assumes that an infinite number of 
repeated measurements on different elements of the ensemble has to be performed to reveal the exact mean value of 
the given observable. In practice only a finite number of measurements can be performed. What is obtained from 
these measurements is a specific set of data indicating number of how many times eigenvalues of given observables have 
appeared (which in the limit of an infinite number of measurements results in the corresponding quantum probability 
distributions). The question is, how to obtain the best a posteriori estimation of the density operator based on the 
measured data. Helstrom Holevo [Q, and Jones j2^ have shown that the answer to this question can be given 
by the Bayesian inference method, providing it is a priori known that the quantum-mechanical state which is going 
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to be reconstructed is prepared in a pure (even though unknown) state. Once this purity condition is fulfiUed, then 
the observer can systematically estimate (i.e. reconstruct) an a posteriori density operator based upon an incomplete 
set of experimental data. This density operator is equal to the mean over all possible pure states weighted by a 
specific probability distribution in an abstract state space with the unique invariant integration measure. It is this 
probability distribution (conditioncd by the assumed Bayesian prior) which characterizes observer's knowledge of the 
System at every moment during the measurement sequence. We note once again that the Bayesian inference has been 
developed for a reconstruction of pure quantum mechanical states and in this sense it corresponds to an averaging 
over a generalized microcanonical ensemble. 

In a real situation one can never design a state-preparation device such that it produces an ensemble of identical 
pure states. What usually happens is that the ensemble consists of a set of pure states, each of which is represented 
in the ensemble with a certain probability (alternatively, we can say that the system under consideration is entangled 
with other quantum-mechanical systems). So now the question is how to use the Bayesian reconstruction scheme 
when the quantum-mechanical system under consideration is in an impure state (i.e., a statistical mixture). To apply 
the Bayesian inference scheme, one has to define exactly three objects: (1) the abstract state space of the measured 
system; (2) the corresponding invariant integration measure of this space; and (3) the prior (i.e., the a priori known 
probability distribution on the given parametric state space). Once these objects are specified one can estimate an a 
posteror density operator after each individual outcome of the measurement has been registered. 

The main purpose of the present paper is to show from an example of the state-reconstruction of a spin-1/2 system 
how the Bayesian scheme of quantum inference developed for a reconstruction of statistical mixtures p4| actually 
works. We will show that this scheme corresponds to a specific averaging over the grand canonical ensemble. Moreover, 
we will show that in the limit of infinite number of measurements the reconstructed density operator is equal to the 
generalized canonical density operator obtained via the Jaynes principle of maximum entropy. In addition, in the case 
the complete observation level (the quòrum of the observables is measured) the generalized canonical density operator 
is equal to the operator obtained via the tomographic measurement. This clearly reveals a logical connection between 
quantum Bayesian inference and quantum tomography. 

The paper is organized as foUows. In Section II we briefly review the Bayesian inference scheme as developed for 
pure states by Jones [ p3[ . In Section III we derive limiting formula for an a posteriori estimated density operator of 
any quantum mechanical system under the prior assumption that the system is in a pure state. Section IV is devoted 
to a description of Bayesian inference for statistical mixtures. Reconstructions of a density operator of a spin-1/2 are 
presented in Section V. We will conclude our paper with general remarks. 



II. BAYESIAN INFERENCE 



The general idea of the Bayesian reconstruction scheme is based on manipulations with probability distributions 
in parametric state spaces. To understand this reconstruction scheme we remind us several definitions and concepts. 
Firstly, it is a space of states of the measured system. The quantum Bayesian method as discussed in the literature 
pï|-p3t is based on the assumption that the reconstructed system is in a pure state described by a state vector |^), 
or equivalently by a pure-state density operator p = |\E')(Í'|. The manifold of all pure states is a continuum which 
we denote as Í2. Secondly, it is the discrete space A of reading states of a measuring apparatus associated with 
the observable O. These states are intrinsically related to the projectors Px- where Ai are the eigenvalues of the 

observable O. 

The Bayesian reconstruction scheme is formulated as a three-step inversion procedure: 

(1) As a result of a measurement a conditional probability 

p(Ó,A,|/5) = Tr(p^^òp). (2-1) 

on the discrete space A is defined. This conditional probability distribution specifies a probability of finding the result 
\i if the measured system is in a particular state p. 

(2) To perform the second step of the inversion procedure one has to specify an a priori distribution po{p) defined on 
the space fl. This distribution describes our initial knowledge concerning the measured system. Using the conditional 
probability distribution p(0, A^jp) and the a priori distribution Po{p) we can define the joint probability distribution 
p{Ó,Xí;p) 

p(Ò,A,;p)=p(Ó,A,|p)po(p), (2.2) 

on the space Í7 ^. We note that if no initial Information about the measured system is known, then the prior po{p) 
has to be assumed to be constant (this assumption is related to the Laplace principle of indifference [p^). 
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(3) The final step of the Bayesian reconstruction is based on the well known Bayes rule p{x\y)p{y) = p{x; y) = 
p{y\x)p{x) , with the help of which we find the conditional probability p{p\0, Xi) on the state space íi: 

Kp|Ó,A.) = -4|^, (2.3) 



from which the reconstructed density operator can be obtained [see Eq.(2.4)]. 

In the case of the repeated 7V-trial measurement, the reconstruction scheme consists of an iterative utilization of the 
three-step procedure as described above. After the A'^-th measurement we use as an input for the prior distribution 
the conditional probability distribution which is an output after the {N — l)st measurement. However, we can 
equivalently define the A^-trial measurement conditional probability p{{ }„|/5) = YÍí^iPÍ^ít ^j\p) ^^'^ applying the 
three-step procedure just once we get the reconstructed density operator 

. InPÜ }N\p)pdn 
JnPil ÏNlPjdn 



where p in the r.h.s. of Eq.(2.4) is a properly parameterized density operator in the state space Q. At this point 
we should mention one essential problem in the Bayesian reconstruction scheme, which is the determination of the 
integration measure dn Q 

The integration measure has to be invariant under unitary transformations in the space íí. This requirement 
uniquely determines the form of the measure. However, this is no longer vàlid when Q is considered to be a space of 
mixed states formed by all convex combinations of elements of the original pure state space íi. Although the Bayesian 
procedure itself does not require any special conditions imposed on the space il, the ambiguity in determination of the 
integration measure has been the main obstacle in generalization of the Bayesian inference scheme for a reconstruction 
of a priori impurc quantum states. 



III. BAYESIAN INFERENCE IN LÍMIT OF INFINITE NUMBER OF MEASUREMENTS 



The explícit evaluation of an a posteriori estimation of the density operator p(A}n) is significantly limited by 
technical difficulties when integration over parametric space is performed [see Eq.(2.4)]. Even for the simplest quan- 
tum Systems and for a relatively small number of measurements, the reconstruction procedure can be technically 
insurmountable problem. 

On the other hand let us assume that the number of measurements of observables Oi approaches infinity (i.e. 
N — + oo). It is clear that in this case mean vàlues of all projectors (P^ ^ ) associated with the observables Oi are 
precis ely know (measured): i.e. 



where = 1. In this case the integral in the right-hand side of Eq.(2.4) can be significantly simplified with the 

help of the foUowing lemma: 

Lemma: Let us define the integral expression 

/(«I, . . . ,q;„_i) = / dxi I dx2..-l dxn-i F{xi, . . . ,Xn-i\ai, . . . ,an-i) (3.2) 

where 



Jo 



Tpl„ I ^ aíN Q2Af cin-iN \a„N rn o\ 

F(a:i,...,a;„_i|Q;i,...,Qí„_i) = — X2 ■■■x^_^ (1 - xi . . . - a;„_i) " (3.3) 
and aj satisfy condition Oj — 1. The integration boundaries yk are given by relations: 



^ Many authors (see, for instance, Ref. identify the prior distribution with the integration measure on the space f2. 

However, the particular form of dn is associated with the topology and the particular parameterization of the space íl rather 
then with some prior Information po(/3) about this system. We will distinguish between these two objects. 
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fe-1 



. . , n — 1 



and B equals to the product of Beta functions B{x,y): 

B = B{an + 1, a„_i + l)S(a„ + a„_i + 1, a„_2 + 2) . . . _B(a„ + a„_i . . . + 1, «i + n - 1). 



ï. The function F{xi, . . . , a;„_i | «i, . . . , q;„_i) in the integral ( |3.2[ ) is a normahzed probabihty distribution 
(n — l)-dimensional volume given by integration boundaries. 
ii. For N — > oo, this probabihty distribution has the foUowing properties: 



{xi) -> a, {xD 



i — 1, 2, 3, . . . , 71 — 1, 



i.e., this probabihty density tends to the product of dcha functions: 



N 



hm F{xx-, ■ ■ .,Xn-i\ai, ■ ■ ■ ,a„_i) = 5{xi - ai)í5(a;2 - «2) • • ■ &{xn-i - an-i)- 



Proof: i. can be derived by a successive apphcation of an equation [see for example Eqs.(3.191)] 

x^-^iu - x)''-^dx ^ ui^+^-^Bin, v). 



ii. can be obtained as a resuh of a straightforward calculations of hmits of certain expressions containin^ 
functions with integer-number arguments. In our calculations we have used the identity 

Í3(n + l,r7i) n 



B{n,m,) n + m 
which is satisfied by Beta functions with integer-number arguments. 



(3.4) 

(3.5) 
in the 

(3.6) 
(3.7) 

(3.8) 
Beta 

(3.9) 



A. Conditional density distribution 

Let us start with the expression for conditional probabihty distribution p({ }jv|/5) for the A^-trial measurement of a 
set of observables Oi. If we assume, that the number of measurements of each observable Oi goes to the infinity then 
we can write: 



p .. p 



(3.10) 



The first product in the right-hand side (r.h.s.) of Eq.(3.10) is associated with each measured observable Oi on a 
given observation level. The second product runs over e igenv alues rii of each observable Oi. 

In what foUows we formally rewrite the r.h.s. of Eq.(3.10): we insert in it a set of (5-function and we perform the 
following integration 



V2 



dxl 



Tt P 



X S 



rii — i 



al.N 



(3.11) 



In Eq.(3.11) we perform an integration over a volume determined by the integration boundaries [see Eq.(3.4)], i.e., 
due to the condition ò P) ^ ^■· tbere is no need to perform integration from —00 to 00. 

At this point we utilize our Lemma. To be specific, firstly we separate in Eq.(3.11) the term, which corresponds 
to the function / given by Eq.(3.2). Then we replace this term by its limit expression (3.7). After a straightforward 
integration over variables x^ we finally obtain an explícit expression for the conditional probabihty p({ }jv^oolp) which 
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we insert into Eq.(^!J) from which we obtain the expression for an a posteriori estimation of the density operator 
p{{ Iat^oo) on the given observation level: 

Pi{ -jfl^U^fl^ [Tr {h,AP) - "}] I P^^- (3-12) 

Here A/" is a normalization constant dctcrmincd by the condition Tr \p(\ In^oo)] ~ 1- 



The mterpretation of Eq.(3.12) is straightforward. The reconstructed density operator is equal to the sum of 



equaUy-weighted pure-state density operators on the manifold Í7, which do sati sfy th e conditions given by Eq.(p.l|) 



[the se con ditions are guaranteed by the presence of 5-functions in the r.h.s. of Eq.( 3.12 )]. In terms of statistical physics 
Eq.( |3.12" ) can be interpreted as an averaging over the generahzed microcanonical ensemble of tho se pu re states which 



do satisfy the conditions on the mean vàlues of the measured observables. Consequently, Eq.( 3.12 ) r epr esents the 
principle of the "maximum entropy" on the generahzed microcanonical ensemble under the constraint (3.1). 

In order to clarify the relationship between the reconstruction procedure based on the Jaynes principio of maximum 
entropy and the quantum Bayesian inference we present in Section V an example of the reconstruction of the state of 
a spin-1/2 systcm. 

IV. BAYESIAN RECONSTRUCTION OF IMPURE STATES 

In classical statistical physics a mixture state is interpreted as a statistical average over an ensemble in which any 
individual realizations is in a pure state. In quantum physics a mixture can be considered a state of a quantum 
system, which can not be completely described in terms of its own Hilbert space, because it is only a part of a more 
complex quantum system. Duc to the lack of information about other parts of this complex system, the description 
of subsystem is possible only in terms of mixtures. 

Let assume quantum system P being entangled with some other quantum system R (reservoir). Let the composed 
system S (composed of P and K) itself is in a pure state \^). The density operator pp of the system P is then 
obtained via tracing over the reservoir degrees of freedom: 

Pp = Ttr{ps\; ps = i*) (í'l. (4.1) 

Once the system is in a pure state, then we can safely apply the Bayesian reconstruction scheme as described in 
Section 2. The reconstruction itself is based only on data associated with measurements performed on the system 
P . When the density operator ps is a posteriori estimated, then by tracing over the reservoir degrees of freedom, we 
obtain the a posterior density operator pp for the system P (with no a priori constraint on the purity of the state of 
the system P). These arguments are equivalent to the "purification" Ansatz as proposed by Uhlmann ]26t . 

To make our reconstruction scheme for impure state consistent we have to chose the reservoir R uniquely. This can 
be done with the help of the Schmidt theorem (see Ref. |^^) from which it foUows that if the composite system S is 
in a pure state jí*) then its state vector can be written in the form: 

M 

|vI/) = ^c.|a,),®|A),, (4.2) 

where \ai)p and are elements from two specific orthonormalized bases associated with the subsystems P and 

R, respectively, and q are appropriatc c omp lex numbers satisfying the normalization condition X^l^íP = 1- The 



maximal index of summation (Af ) in Eq.(4.2) is given by the dimensionality of the Hilbert space of the system P. In 
other words, when we apply the Bayesian method to the case of impure states of M-level system, it is sufficient to 
"couple" this system to an M-dimesional "reservoir" . Due to the fact that we measure only observables of the first 
subsystem P particular form of states lA)^ of the second subsystem R does not affect results of the reconstruction. 

V. SPIN-1/2 RECONSTRUCTION 

We assume an ensemble of spins-1/2 in an unknown state described by the density operator 

p{eA)^l{i + r^) (5.1) 
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where f = {sm9cos(j),sm9s'm(j>,cos0); (f> G (0,27r), 6 G (0, tt), and ï is the unity operator. The Pauli spin operators 
à in the matrix representation in the basis |0}, |1) of the eigenvectors of the operator az do read 



01\ ^ _ f -i\ . _/10 

1 j ' '^y^[iol· ^ l -1 



(5.2) 



To determine completely the unknown state, one has to measure three hnearly independent (e.g., orthogonal) 
projections of the spin-1/2. One possible choice of the complete set of observables (i.e., the quòrum associated 
with the spin-1/2 arc spin projections for three orthogonal directions represented by the Hermitian operators: 

= y, i = x,y,z (5.3) 
In what foUows we will consider three observation levels defined as Oa = {sz}, Ob = {sz,Sx} and Oc — 

{^Zí — ^comp- 

ili order to apply the Bayesian inference scheme for impure states, we have to "purify" ||2^ the quantum-mechanical 
system as discussed in Section IV. In the particular case of the spin-1/2 it means that we have to consider a system 
of two spins-1/2 on of which does play the role of a reservoir. In this case the corresponding two-spins-1/2 density 
operator can be parameterized as 



10 1 f<i)CT (g) f(2)(T 



p{<^,ÍJ,(t>i,di,h,d2) = —. \ -. hcosa H (5.4) 



f(l)CT(X)l l®f(2)(T 



sin a cos ip 



4 4 
where i/', </'ij ^2 G (0, 27r); a, 61,62 £ (0, tt) and 



sin a sin ■0 



fc<j' = (sin 0j , — cos , 0) ; 

íü' = (cos 6j cos (f)j , cos 6j sin ipj , — sin dj); 

f'^' = (sin 9j cos (pj , sin 6j sin (pj , cos 6j ) . 



(5.5) 



Once we have parameterized the state space Í7 we have to find the invariant integration measure dçi . We have derived 
this measure earlier ÍE4| and it reads: 



dn = cos^ a sin a sin 01 sm62dad^pd(f>iddid<j)2dd2. 



(5.6) 



A set of projectors associated with the observables a ^ and a ^ [in what foUows we use the notation such that the 

(5.7) 



position of the operator to the left (right) of the symbol (E) is associated with the first (second) spin-1/2] 

(i + SCT,) 



{1 + sài) ^ 



p 



1 



- 1 (g) 1 ài®àj 

.(1) .(2) = — h S — ^. 



The corresponding conditional probabilities do read 
p(s,sr'|p(a...)) = l + .^rr'; 



p(.,5f' |p(a...)) = Í-t-/°'(")^<^ 



2 ■'^^ ' 



(5.8) 



p(s, èf èf |p(a ...)) = ^ + + s^-^ [(fc^' kf - l^P If ) cos 7^ - (fcf' If + fcf ) sin ^] • 

We remind ourselves that we do consider only measurements performed on the first spin described by the observables 
(t'^' . After the Bayesian reconstruction of the composed system is performed then the "reservoir" degrees of freedom 
are traced out. The resulting density operator describes an o posteriori estimation of the density operator of a 
two-level system with no reference on the a priori assumption about the purity of the spin. 

Instead of analysing estimated density operators after a finite number of measurements performed over the spin-1/2 
(i.e. measurements performed over a finite number of elements of the ensemble) we focus our attention on results of 
the reconstruction in the limit of infinite number of measurements. 
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A. Observation level Oa = {sz} 



On the observation level O a only the spin component Sz is nieasured. This kind of the measurement can be 
performed with the help of one Stern-Gerlach apparatus. With the help of the data obtained in a large (infinite) 
number of measurements we can express the density operator of the spin-1/2 under consideration as (the trace over 
the reservoir degrees of freedom has already been performed) 

P=-Jj\ y'dy^ sin^id^i (5((<T<^') -ycosí?i)(Í+2/cos0ià,), (5.9) 

where the variable a is substituted by y = cosa. When we perform integration over the variable y we obtain the 
expression 

P=\f\ n | (í + {^'^à.), (5.10) 

J\ J^i cos^ t^i I cos t^l I 

where the integration is performed over the region 

C :={(), tt} suchthat | cos6li| > |((T<'')|. (5.11) 
After we perform the integration over 9i we obtain for the a posteriori estimation of the density operator the expression 

P=l{'^ + {àz)àz). (5.12) 
which is identical to the one obtained via the Jaynes principle of maximum entropy [p8| . 

B. Observation level Ob = {Sz,Sx} 

Let us extend the observation level Oa and let us assume the measurement of two spin projections Sz and Sx (i.e. 
two Stern-Gerlach apparatuses with fixed orientations are employed). In the limit of infinite number of measurements 
one can express the Bayesian estimation of the density operator of the spin-1/2 on the given observation level as (here 
the trace over the "reservoir" spin has already been performed): 

P^JfJ y^^yj Sïn9id9iJ d(j)i í((<tÍ'' ) - y cos 6'i)(5(((tÍ'' ) - y sin 6*1 cos 0i) 

X (1 + j/siné*! cos^iàa; + ysm9i sin^iCT^ -|- ycos^itr^). (5.13) 
When we perform integration over the variable y we find 

1 í^^ f ,^ siné»! 



p=- dçbi d9i . J((à<^')-tangicos0i(ài^')) 

M Jo Jc cos^ 9i I cos 9i \ 

X (i + ((TÍ'')tané'iCOS(/)i<Ta; -|- ((tÍ'') tané»! sin(/)i(Ty -|- (<tÍ'' )ò·z). (5.14) 



The integration over the variable 0i in the right-hand side of Eq. ( p.l4| ) results into the foUowing expression 



2 

uf ^^lè 7-^(1 + {àí'')àx + {àí'')tan9ism4'^ày + {àí'')àz), (5.15) 

^ JC" cos^ 9i I sin (pi ' I 



where the integration is performed over the region 

£":={0,7r} suchthat \cos9i\ > \{à'z^)l and | tan é'l | > ^-^^ 

{a'z') 



(5.16) 
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The sum in Eq.(5.15) is performed over two vàlues (j)"-^^ oi the variable (pi which are equal to two solutions of the 
equation 



cos <pi 



(TÍ'*)tan6li 



(5.17) 



Due to the fact that the term in front of the operator ct^^' is an odd function of (f)'"/^ , we can straightforwardly perform 
integration over 9i and we find the expression for the reconstructed density operator 



(5.18) 



which again is exactly the same as if we perform the reconstruction with the help of the Jaynes principle 



C. Complete observation level Oc ~ {sí''\ s^x \ s'jj'^} 

Let us assume now that the measurement is performed with the help of three Stern-Gerlach apparatuses, each of 
which are measuring three spin component s^·' (i = x,y,z). On this complete observation level the expression for 
the Bayesian estimation of the density operator of the spin- 1/2 in the limit of infinite number of measurements can 
be expressed as (here again we have already traced over the "reservoir" degrees of freedoms) : 

P=JJ j y^dyj svLíOidOij (i0i(5(((7Í''} - y cos 6'i)(5(((74''} - y sinó'l cos (/)i) 



X ò{{ày ^) — y sinó'l sin(/)i)(ï + ysiné'i cos^ià^: + y sinó'l sin0i(7j^ + ycosOiàz 



(5.19) 



We can rewrite Eq.( 5.19 ) as 



2 

Tr í ^^lE ,J ^ ((4^')-tangiSÍn^<^·)(ài^')) 
■1^ Jc" . 1 cos^ é'l! smç!)^-'^! ^ ^ 



i + {àí'')à^ + (4^') tan 01 smc^Y^ày + {ài'^)à 



(5.20) 



are defined by Eqs.( 5.16 ) and ( 5.17 ), respectively. Now the integration over parameter 9i can be 



where £" and 

easily performed and for the density operator of the given spin-1/2 system we find 



i(i + (<7<^'}a, + (4^')à, + (à<^>)à,). 



(5.21) 



The density operator (5.21) obtained with the help of Bayesian inference scheme is equal to that one which foUows from 
the Jaynes principle of the maximum cntropy. Moreover this same result can be obtained from other reconstruction 
schemes, such as the discretc quantum tomography |^] (see also and p9|). 

The simple example of the spin-1/2 reconstruction reveals deep conceptual relationship between the quantum 
Bayesian inference and the Jaynes principle of the maximum entropy. To understand this relationship more clearly we 
turn our attention once again to the a priori assumption under which we have performed the reconstruction. We have 
assumed that the spin-1/2 can be in an impure state (therefore we have applied the "purification" procedure). As a 
consequence of this assumption the von Neumann entropy of the reconstructed density operator may be larger than 
zero, which is equivalent to the fact that the mean vàlues of the observables òi {i — x,y, z) do fulfiU the condition 



i^xf + {àyf + {àzf < 1, 



(5.22) 



i.e. the reconstructed state can be expressed as a point either on or inside the Poincare sphere. 

On the contrary, if it is a priori assumed that the reconstructed state is a pure one (see, for instance, works by Jones 
p3|), then the Bayesian reconstruction in the limit of infinite number of measurements on the complete observation 
level results in the reconstructed density operator which can be expressed as 
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27r 



P — I sm9d0 S{{àz) — cos 9)6 {{àx) — sin 9 cos (jj) 6 {{ày) 







X (i + sin 9 cos (j>(Tx + sin 9 sin (p&y + cos 9àz)- (5.23) 

The integral in the right-hand sidc of the equation can only be performed if the mean vàlues of the observables under 
consideration do fulfill the purity condition 

{àxf + {àyf + {cTzf = 1. (5.24) 



Providing the purity condition holds then from Eq.(5.23) we obtain for the density operator the expression ( ^.21 ), 
otherwise the reconstruction scheme fails. 

The limit formulae for the Bayes inference have an appealing geometrical interpretation. For example the three 



(5-functions in Eq.(5.23) correspond to three specific orbits on the Poincare sphere, each of which is associated with a 
set of pure states which posses the measured value of a given observable s^. This corresponds to an averaging over the 
generalized microcanonical ensemble of pure states having the measured mean vàlues of the three observables. The 
reconstructed density operator then describes a point on the Poincare sphere which coincides with an intersection of 
these three "orbits". Consequently, if the three orbits have no intersection, the reconstruction scheme fails, because 

there does not exist a pure state with the given mean vàlues of the measured observables. 

On the other hand, the form of the a posteriori expression for the reconstructed density operator ( 5.19| ) obtained 



with the help of the Bayesian inference with no a priori restriction on the purity of the reconstructed state reveals 
that the estimated density operator can be obtained as a result of averaging over all points on and inside the Poincare 
sphere. This averaging over the generalized canonical ensemble is equal to the maximization of the vonNeumann 
entropy as assumed by Jaynes. 



VI. CONCLUSIONS 



In the paper we have analyzed in detall the logical connection between three different reconstruction schemes: 
(1) If measurements over a finite number of elements of the ensemble are performed then one can obtain the a 
posteriori estimation of the density operator with the help of the Bayesian inference. If nothing is know about the 
reconstructed state one has to assume a constant prior probability distribution on the parametric state space under 
the assumption that the system is in a statistical mixture. (2) As soon as number of measurements becomes large 
then the Bayesian inference scheme becomes equal to the reconstruction scheme based on the Jaynes principle of the 
maximum entropy, i.e., in the limit of infinite number of measurements a posteriori estimated density operator fulfills 
the condition of the maximum entropy. Consequently, it is equal to the generalized canonical density operator. (3) If 
the quòrum of observables is measured, then the generalized canonical operator is equal to the "true" density operator 
of the system itself, i.e. the complete reconstruction via the MaxEnt principle is performed. It is the question of 
technical convenience which reconstruction scheme on the complete observation level is utilized (for instance, quantum 
tomography can be used), but the fact is that all of them can be formulated as a maximization of the entropy under 
given constraints. 
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